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A broad absorption feature in the X-ray spectrum of the 
isolated neutron star RBS1223 (1RXSJ1 30848.6+21 2708) 
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Abstract. X-ray spectra of the isolated neutron star RBS1223 obtained with the instruments on board XMM-Newton in 
December 2001 and January 2003 show deviations from a Planckian energy distribution at energies below 500 eV. The spectra 
are well fit when a broad, Gaussian-shaped absorption line with a = 100 eV and centered at an energy of 300 eV is added to 
an absorbed blackbody model. The resulting equivalent width of the line is — 150 eV. However, the spectral resolution at these 
low energies of the EPIC detectors and the lower statistical quality and restricted energy band of the RGS instruments are not 
sufficient to exclude even broader lines at energies down to 100 eV or several unresolved lines. The most likely interpretation 
of the absorption feature is a cyclotron absorption line produced by protons in the magnetic field of the neutron star. In this 
picture line energies of 100—300 eV yield a magnetic field strength of 2— 6x 10 13 G for a neutron star with canonical mass and 
radius. Folding light curves from different energy bands at a period of 10.31 s, which implies a double peaked pulse profile, 
shows different hardness ratios for the two peaks. This confirms that the true spin period of RBS1223 is twice as long as 
originally thought and suggests variations in cyclotron absorption with puls e phase. We a l so pro pose that changes in photo- 
electric absorption seen in phase resolved spectra of RX J0720.4— 3125 by ICropper et alj 1200 lh . when formally fit with an 
absorbed blackbody model, are caused instead by cyclotron absorption varying with pulse phase. 
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1. Introduction 

The soft X-ray source 1RXS J130848.6+212708 was discov- 
ered in the ROSAT all-sky survey data and further stud- 
ied as part of the ROSAT Bright Survey program (RBS, 
source number 1223) to optically identify the brightest ^2000 
high-galactic latitude sources (Schwope et al. 2000). Based on 
similarities to the two best known cases RX J 1856.4—3754 
and RX J0720.4— 3125 and the lack of an optica l counter- 
part brighter than B^26 m , ISchwope et alJ {l999) proposed 
1RXS J130848.6+212708 = RBS 122 3 as isolated neutron star 
(INS). Using the Chandra observatorv lHambarvan et all J2002 ) 
observed RBS 1223 on June 24, 2000 with the ACIS instrument 
and discovered pulsations in the X-ray flux, initially indicat- 
ing a neutron star spin period of 5.16 s. A probable optical 
counterpart in the 9 0% Chandra error circle was reported by 
Kapl an et al . (2002a) using a very deep observation from the 
Hubble Space Telescope. The optical brightness of msoccD 
= 28.6'" yields an X-ray to optical flux ratio of log(f x /f opt ) = 
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4.9. The optical flux is a factor of about 5 above the extrapo- 
lation of the blackbody fit to the X- ray spectrum, similar as 
observed for RX J 1 856.4— 3754 ("e. e. iBurwitz et alJl20ol and 
RX J0720.4— 3125 JPavlov et all2002l) . 

This confirms RBS 1223 as member of a group of X-ray 
dim INSs which share similar properties (for a recent review 
see lHaberll 12003): soft blackbody-like X-ray spectrum, radio 
quiet and no known association with a supernova remnant. 
Four objects are X-ray pulsars with neutron star spin periods 
in the range 8.39 to 22.7 s (the XMM-Newton observations of 
RBS 1223 revealed a peak in the power spectrum at half the fre- 
quency of the main peak which suggests that the genuine spi n 
period is 10.31 s with a double peaked profile. lHaberll E"003). 
All available spectra obtained by the ROSAT PSPC were con- 
sistent with Planckian energy distributions with blackbody 
temperatures kT in the range 40 — 100 eV and little attenua- 
tion by interstellar absorption. More recent obse rvations of the 
two b rightest objects RX J18 56.4-3754 dBurwitz et alJl200ll 
2003) and RX .10720.4-3 125 dPaerels et alJ200lUPavlov et alJ 
2002) were performed using the low energy transmission grat- 
ing (LETG) aboard Chandra and the reflection grating spec- 
trometers (RGS) of XMM-Newton. Also at high spectral reso- 
lution both sources show featureless spectra which can best be 
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Table 1. XMM-Newton observations of RBS1223. 



Time(UT) XMM Read-out Count rate Exp. 

Start End Detector Mode [s -1 ] [ks] 

2001 Dec. 31, Satellite Revolution 377 (AOl) 

03:08 08:41 MOS1/2 Full Frame 0.49/0.52 19.8 

03:32 08:42 pn Small Window 2.47 19.0 

03:02 08:42 RGS1/2 Spectro+Q 0.09/0.07 20.4 

2003 Jan. 1, Satellite Revolution 561 (CAL) 

06:13 14:14 MOS1/2 LargeWindow 0.47/0.48 28.7 

06:36 14:14 pn Full Frame 2.39 27.0 

06:12 14:15 RGS1/2 Spectro+Q 0.10/0.08 28.9 

Mean count rates are given for the energy band used for the spectral 
analysis (EPIC: 0.13-1.3 keV, RGS: 0.32-0.90 keV). The thin filter 
was used in all EPIC cameras during both observations. 



modeled by a Planckian spectrum with kT of ^63 eV and ^86 
eV, respectively. 

The pulsars among the X-ray dim INSs may allow an es- 
timate of the strength of their magnetic field if their long- 
term spin period increase can be interpre ted in terms of dipolar 
losses. For the pulsar RXJ0720.4-3125 lkaplan et alJ (|2002b) 
derived an upper limit for the period derivative of 3.6xl0 -13 
ss" 1 from the analysis of ROSAT, S AX and Chandra d ata and 
including first XMM-Newton data, IZane et all (2002) deter- 
mined possible values between 3—6 x 10~ 14 ss" 1 . This implies 
that the magnetic field strength of RX J0720.4-3125 could be 
as high as a few 10 13 G. For such magnetic field strengths an 
absorption feature at the proton cyclotron resonance energy, 
E cp = 63g r (B/10 13 G) eV (g r = (l-2GM/c 2 R) 1 / 2 with M and 
R the neutr on star mass and radius, respectively), ma y become 
detectable dZane et alJl200ll IZavlin & Pavlovll2002l) with de- 
tectors sensitive at soft X-rays down to 100 eV. 

In this letter we present the results from a spectral analysis 
of two XMM-Newton observations of RBS1223. The spectra 
show deviations from a Planckian energy distribution which 
can best be modeled by a broad absorption feature at low en- 
ergies. We discuss this feature as most likely being caused by 
proton cyclotron resonance absorption. 



2. XMM-Newton observations 

XMM-Newton Jjansen etall200ll) observed RBS1223 on two 
occasions; in the guest observer program (AOl, PI: Schwope) 
and for calibration purposes at the beginning of the year 
2003. Here we utilize the data collected with the European 
Photon Imaging Cameras (EPIC s) and the Refl ection Grating 
Spectrometers (RGSs, Iden Herder et alJ boOlh for a spectral 
analysis. Two EPIC cameras are based on MOS (EPIC-MOS1 
and -MOS2 iTurner et alJ l200ll) and one on pn (EPIC-pn, 
Struderetal .2001) CCD detectors and are mounted behind the 
three X-ray multi-mirror systems. The details of the XMM ob- 
servations are summarized in Tabled 



3. X-ray spectra 

The data from the EPIC instruments were processed using 
SAS5.4.1 to produce event files for pn, MOS1 and MOS2. 
RGS spectra and response files were created using rgsproc of 
SAS5.4.1 with standard settings. The EPIC spectra were ex- 
tracted from circular regions with 30" radius around the tar- 
get position and nearby background regions using single-pixel 
events (pattern 0) from pn and pattern 0—12 events from MOS 
data. Detector response matrix files were used as available 
from the EPIC-pn calibration group in February 2003 (ver- 
sion 6.5 which will be implemented in the next SAS release 
after SAS5.4.1) and were created for MOS using the arfgen 
and rmfgen tasks of SAS5.4.1. MOS spectra from recent ob- 
servations suffer from a considerable change in the low en- 
ergy re-distribution which is not calibrated yet. Therefore, we 
use MOS spectra from the first observation in Dec. 2001, only. 
The pn spectra from the two observations, obtained in different 
read-out modes, are in agreement within the current calibra- 
tion uncertainties and also are consistent with the MOS spectra 
from the Dec. 2001 observation. 

The spectra were simultaneously fitted with various mod- 
els, all including photo-elect ric ab s orption with element abun- 
dances from Anders & Grevesse ( 1989) and cros s-sections 
from Balucinska-Church & McCammon i 19921) and lYanetalJ 
dl998l) using XSPEC version 11.2.0. To account for cross- 
calibration uncertainties between the instruments and the fact 
that the pn spectra were not corrected for point spread function 
losses (the SAS5.4.1 task arfgen is not compatible to the re- 
sponse files we used here) the relative normalization between 
the instruments/observations were allowed to vary individually. 
The pn fluxes were corrected manually by a factor 1 .2 to correct 
for the PSF loss according to the XMM-Newton users hand- 
book. In addition uncertainties in the low-energy re-distribution 
were accounted for by allowing individual column density val- 
ues for pn, MOS 1 and MOS2. The RGS spectra show relatively 
large differences in continuum parameters which were then fit 
for each spectrum separately. 

As first model an absorbed blackbody was used, however, 
the combined fit to the XMM spectra of RBS1223 is not ac- 
ceptable with Xr °f 5-2 for 639 degrees of freedom (dof). The 
spectra together with the fit are shown in Fig. [2 and the fit pa- 
rameters summarized in Table|2] Because of the bad fit we do 
not give errors and fluxes for this model. In particular, the EPIC 
spectra with best statistical quality yield high \ 2 values which 
we report for each instrument separately in Table [2] Since the 
combined fit is dominated by the high statistics of the EPIC-pn 
spectra we allowed all model parameters for the RGS spectra 
to be fit freely, i.e. the RGS fit is completely independent from 
the pn fit. This allows to judge if the RGS spectra on their own 
can be well fit by a blackbody model. 

We a lso tried non-magnetic neutron star atmosph ere mod- 
els (e.g. Gans icke et alJ2 002: Zavlin & Pavlov 2002). Iron and 
solar mixture atmospheres cause too many absorption features 
and deviations from a blackbody model in particular at ener- 
gies between 0.5 and 1 .0 ke V which are not seen in the mea- 
sured spectra. On the other hand the spectrum of a pure hy- 
drogen model is similar in shape to that of a blackbody and 



Table 2. Spectral fit results. 
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Model A: phabs*(bbody) 




Model B: phabs*(bbody+gaussian) 






Obs. -Instrument 


kT 


N H 


X 2 /dof 


kT 


N H 


Elinc 


EQW 


X 2 /dof 


Flux (1) 






[eV] 


[10 20 cm" 2 ] 


per instr. 


[eV] 


[10 20 cm" 2 ] 


[eV] 


[eV] 


per instr. 


[erg cm" 2 


r 1 ] 


AOl-pn 


95.1 


7.1 


2138/309 


85.8±0.5 


4.1±0.1 


290±5 


-148 


589/307 


3.65x10" 


12 


CAL-pn 


= 1(2) 


= 1 




=1 


=1 


=1 


-148 




3.44x10" 


12 


AOl-MOSl 


= 1 


4.9 


845/109 


=1 


2.9±0.2 


302±2 


-159 


156/108 


3.33x10" 


12 


A01-MOS2 


= 1 


5.3 




=1 


3.2±0.2 


=3 


-159 




3.62x10" 


12 


AOl-RGSl 


103 


7.1 


349/221 


82.2±2.4 


5.0±0.6 


= 1 




318/221 


2.62x10" 


12 


A01-RGS2 


109 


=5 




85.2±2.6 


=5 


= 1 






2.50x10" 


12 


CAL-RGS1 


112 


=5 




87.8±1.6 


=5 


= 1 






2.66x10" 


12 


CAL-RGS2 


108 


=5 




87.6±1.5 


=5 


= 1 






2.62x10" 


12 


ACIS 


81.1 


19 


253/38 


87.8±1.0 


<1.6 


290fixed 




91/38 


3.34x10" 


12 



Observed flux 0.1—2.4 keV; ^ "=n" denotes fit parameter is linked with parameter in line n 




0.1 2 0.5 1 0.1 0.3 0.5 1 

Channel Energy (keV) Channel Energy (keV) 



Fig. 1. Blackbody model fits to EPIC-pn (upper pair), EPIC-MOS (middle pair) and RGS spectra of RBS1223. The four RGS 
spectra were combined in the plot for clarity. While the pure blackbody model fit (left) is unacceptable, including a broad 
Gaussian absorption line at ^300 eV (right) can reproduce the data. The residuals (bottom panels) show consistent behavior for 
all instruments. 



does not fit the data either. It also results in a much lower effec- 
tive temperature of abo ut 30 eV w hich would predict a far too 
high optical flux (seelPavlov et alJfT996l) . Using the models of 
lGansickeetaT1ll2002l) xl values of 12, 29 and 12 (for the com- 
bined EPIC spectra) are obtained for Fe, solar mixture and H 
atmospheres, respectively. A two-temperature blackbody fit re- 
sults in Xr = 6-6 with kTi = 33 eV and kT2 = 94 eV. The strong 
residuals at energies below 500 eV remain. Such residuals are 
not seen in fits to the spectra of RX J 1856.4—3754 with typi- 
cal reduced Xr values of 1.5 using the same version of spectral 
response files, clearly excluding calibration problems as origin. 

An acceptable fit was found when including a broad ab- 
sorption line in the model. Because of the low line energy and 
the large line width the low resolution of the EPIC detectors 
and the limited band pass of the RGS instruments do not allow 
to independently vary energy and width of the line. A broad 
line centered at the lowest energies reachable by the detectors 
(~120 eV) leaks only partially into the sensitive energy band 
and can therefore not be distinguished from a somewhat nar- 
rower line at somewhat higher energy. All we can say is that 



the absorption feature is extended up to ^500 eV. Models with 
an absorption line at 300 eV with a = 100 eV or a line at 200 eV 
with a = 150 eV can not be distinguished. Also multiple lines 
at e.g. 150 eV and 300 eV which are correspondingly narrower 
can not be exclud ed. However, too broad lines are probably un- 
realistic (see e.g. lZane et al 1 l200ll if interpreted as cyclotron 
resonance) and we therefore use in our fits a single absorption 
line of Gaussian shape with fixed a = 100 eV. This yields a best 
fit with Xr = 1-67 for 636 dof (Fig.[0 with parameters listed in 
Table 13 The line energy was allowed to be fit individually for 
the two EPIC instruments but is consistent within the errors. 
Because the sensitive energy band of the RGS instruments cov- 
ers the broad line only partially, we forced all line parameters 
to be the same in the fit for the pn and all RGS spectra simulta- 
neously. For the same reason we do not give values for the line 
equivalent width (EQW) for the RGS spectra. The x 2 obtained 
from the RGS spectra for this model is reduced by 3 1 (without 
adding additional dofs) as compared to the blackbody model 
which may be accepted for RGS if treated by its own. 
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Fig. 2. Blackbody model fits with (bottom) and without absorp- 
tion line to the ACIS-S spectrum of RBS1223. 

We r e-analyzed the ACIS-S s pectrum of RBS1223, pub- 
lished bv lHambarvan et all ll20Q2ll . using the re-processed data 
available in the Chandra archive and new calibration which al- 
lows a spectral analysis down to 350 eV. We fit blackbody with 
and without absorption line to the spectrum, with all line pa- 
rameters fixed at the values derived from the fit to the XMM- 
Newton spectra. The fits are shown in Fig. |2] and the model pa- 
rameters reported in Table|2] The results are in good agreement 
to those obtained from the XMM-Newton data. 

4. X-ray pulsations 

To investigate spectral changes with pulse phase, light curves in 
the energy bands 0.12—0.5 and 0.5 — 1.0 keV were folded with 
a pulse period of 10.31 s. This long period with double peaked 
profile is suggeste d by the pr esence of a second peak in the 
power spectrum faaberllE o03). The different relative strength 
of the two peaks in the two energy bands and the related sig- 
nificant difference in hardness ratio (Fig. [3} strongly supports 
a 10.31 s spin period of RBS1223, twice as long as originally 
thought. 

5. Discussion 

Spectral analysis of the INS RBS1223 using data obtained by 
the X-ray instruments aboard XMM-Newton revealed signifi- 




0.5 1 1.5 



Phase (10.312556 s) 

Fig. 3. Pulse profile of RBS1223 in the 0.12-0.5 keV (soft) 
and 0.5 — 1.0 keV (hard) energy bands, together with the ratio 
hard/soft, obtained from the EPIC-pn data of the Jan. 2003 ob- 
servation. 

cant deviations from a Planckian energy distribution at energies 
below ^500 eV. Also, non-magnetic atmosphere models yield 
unacceptable fits. We find that a model including a broad ab- 
sorption line on top of the blackbody continuum can reproduce 
the observed spectra. For a Gaussian line with fixed er= 100 eV 
we derive an energy of the line of ^300 eV and an equivalent 
width of — 150 eV. However, the spectral resolution of the EPIC 
instruments at these energies and the restricted energy coverage 
of the RGS are not sufficient to exclude line energies as low as 
100 eV or multiple lines which are unresolved. 

Cyclotron resonance absorption features in the 0.1 — 1.0 
keV band are expected in spectra from magnetized neutron 
stars with field strengths in the range of 10 10 — 10 11 G or 
2x 10 13 — 2 x 10 1 4 G if caused by e lectrons or protons, res pec- 
tively (see e.g. IZane et all 1200 it Izavlin & Pavlovl l2002l) . A 
strong ma gnetic field for RBS122 3 of ~2xl0 14 G was first 
derived by Hamba rvan et alJ d2002t) from a first estimate of the 
spin period derivative based on a ROSAT and a Chandra ob- 
servation. The first of the two XMM-Newton observations de- 
scribed her e, however, p ut this into question (for preliminary 
results see lHaberll E003> and placed only an upper limit of 
10 14 G. A detailed temporal analysis of the two XMM-Newton 
observations and a re-analysis of the Chandra data will be pub- 
lished in an accompanying paper (Hambaryan et al. in prepara- 
tion). 

Recently, several radio pulsars with magnetic field 
strength of a few 10 13 G with similar long spin pe- 
riod as RXJ0720.4-3125 and RBS12 23 were discovered 
JCamilo et all l200fj iMorris et all l2002h . If one assumes a 
comparable field strength for RBS1223 which may be jus- 
tified if the ROSAT discovered radio-quiet IN Ss are cases 
where the radio beam doe s not cross the Earth jKaplan et alJ 
2002b; iMotch et alJEool . but are otherwise similar objects, 
one would favor cyclotron resonance absorption due to pro- 
tons to be responsible for the features seen in the spectra of 
RBS1223. Then our measured line energies between 100 and 
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300 eV translate into magnetic field strengths of 2—6 x 10 G 
for a neutron star with canonical mass (1.4 Mq) and radius (10 
km). 

By folding the EPIC-pn data of RBS1223 in two energy 
bands with a period of 10.31 s we see hardness ratio changes 
with pulse phase. The two pulses show different hardness 
which may be explained by variations in cyclotron absorption 
between two magnetic poles. Detailed phase resolved spec- 
troscopy which exceeds the scope of this letter, is in progress 
and will be reported elsewhere. 

Variations in the hardness ratio with pulse phase were also 
reported from RX J0720.4— 3125. Formal fits using a black- 
body spect rum attenuated by p hoto-electric absorption per- 
formed by Cropp er et alJ (|2001 ) to the pulse-phase resolved 
EPIC-pn spectra around hardness maximum and minimum 
yielded changes in column density, but little variations in black- 
body temperature. This could be explained in terms of energy- 
dependent beaming effects with softer photons more strongly 
beamed than harder photons. From their polar cap modeling of 
the spin pulse profile the authors conclude that the derived ge- 
ometries are at best marginally compatible with surface temper- 
ature variation but more satisfactory explained by an accretion 
model. However, the possibility of accretion of interstellar mat- 
ter is now most likely exclud ed by the measure d high proper 
motion of RX J0720.4-3125 Jlvlotch et all2003l) . On the other 
hand, the variations with pulse phase could be naturally ex- 
plained by cyclotron absorption, because we may view surface 
areas with different magnetic field strength and orientation dur- 
ing the rotation of the neutron star. A cyclotron absorption line 
at the low energy end of the sensitive energy band of the EPIC 
instruments would cause deviations from a Planckian shape at 
these energies which, so far, were not recognized due to lack of 
statistics and/or energy resolution and preliminary calibration 
of the new generation instruments. 

IZane et all J200ll) calculate spectra emerging from ultra- 
magnetized neutron stars covered by non-degenerate, pure hy- 
drogen plasma. They explored ranges in magnetic field strength 
from 10 13 G to 10 15 G and in luminosity from 10 34 erg s _1 
to 10 36 erg s _1 . Although the luminosity of RBS1223 may 
only be around lfj 31 ~ 32 erg s _1 their results for the line EQW 
shows little dependence on luminosity and may be compared 
to the value of 150 eV we find for RBS1223. Their predictions 
for cyclotron line widths and equivalent width and the over- 
all blackbody-like continuum for the case of a dipolar mag- 
netic field (which broadens the line due to B field variations 
in magnitude and direction along the neutron star surface) are 
consistent with our results for RBS1223. This suggests, that 
at least for RBS1223 and probably also RX J0720.4— 3125 a 
H atmosphere exists on these neutron stars and argues against 
a condensed matter surface as it was proposed to explain the 
combined spectral energy d istributio n in the o ptical and X- 
ray bands jBurwitz et al 1200*31 IZane et alj f2003). However, an 
optically thick H atmosphere produces too much optical flux 
(Pav lov et alJ ll996). a problem which may be reduc e d by t hin 
atmosphere models as was applied by iMotch et alJ J2003I) to 
RXJ0720.4-3125. 
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